Genomic epidemiology reveals multidrug resistant plasmid spread between Vibrio cholerae lineages in Yemen.

Since 2016, Yemen has been experiencing the largest cholera outbreak in modern history. Multidrug resistance (MDR) emerged among Vibrio cholerae isolates from cholera patients in 2018. Here, to characterize circulating genotypes, we analysed 260 isolates sampled in Yemen between 2018 and 2019. Eighty-four percent of V. cholerae isolates were serogroup O1 belonging to the seventh pandemic El Tor (7PET) lineage, sub-lineage T13, whereas 16% were non-toxigenic, from divergent non-7PET lineages. Treatment of severe cholera with macrolides between 2016 and 2019 coincided with the emergence and dominance of T13 subclones carrying an incompatibility type C (IncC) plasmid harbouring an MDR pseudo-compound transposon. MDR plasmid detection also in endemic non-7PET V. cholerae lineages suggested genetic exchange with 7PET epidemic strains. Stable co-occurrence of the IncC plasmid with the SXT family of integrative and conjugative element in the 7PET background has major implications for cholera control, highlighting the importance of genomic epidemiological surveillance to limit MDR spread.

Genomic epidemiology of Vibrio cholerae during a mass vaccination campaign of displaced communities in Bangladesh.

Ongoing diarrheal disease surveillance throughout Bangladesh over the last decade has revealed seasonal localised cholera outbreaks in Cox's Bazar, where both Bangladeshi Nationals and Forcibly Displaced Myanmar Nationals (FDMNs) reside in densely populated settlements. FDMNs were recently targeted for the largest cholera vaccination campaign in decades. We aimed to infer the epidemic risk of circulating Vibrio cholerae strains by determining if isolates linked to the ongoing global cholera pandemic ("7PET" lineage) were responsible for outbreaks in Cox's Bazar. We found two sublineages of 7PET in this setting during the study period; one with global distribution, and a second lineage restricted to Asia and the Middle East. These subclades were associated with different disease patterns that could be partially explained by genomic differences. Here we show that as the pandemic V. cholerae lineage circulates in this vulnerable population, without a vaccine intervention, the risk of an epidemic was very high.

Phylogenomics reveals insights into the functional evolution of the genus Agrobacterium and enables the description of Agrobacterium divergens sp. nov.

The genus Agrobacterium was initially described as mainly phytopathogenic strains. Nowadays, the genus includes phytopathogenic and non-phytopathogenic bacteria that are distinctive among the Rhizobiaceae family. Recently we have isolated two closely related strains, LMG 31531T and LMG 31532, from soil and plant roots, respectively. Both strains differ from previously reported species based on the genomic and phenotypic data. A. arsenijevicii KFB 330T and A. fabacearum LMG 31642T showed the highest 16S rRNA similarity (98.9 %), followed by A. nepotum LMG 26435T (98.7 %). A clear genomic feature that distinguishes LMG 31531T and LMG 31532 from other Agrobacterium species is the absence of a linear chromid. Nevertheless, typical values of the core-proteome Average Amino Acid Identity (cpAAI > 85 %) and 16S rRNA gene sequence similarity (>96 %) when compared to other members of the genus confirm the position of these two strains as part of the Agrobacterium genus. They are therefore described as Agrobacterium divergens sp. nov. Besides, our comparative genomic study and survey for clade-specific markers resulted in the discovery of conserved proteins that provide insights into the functional evolution of this genus.

Genomic analysis provides novel insights into diversification and taxonomy of Allorhizobium vitis (i.e. Agrobacterium vitis).

BACKGROUND: Allorhizobium vitis (formerly named Agrobacterium vitis or Agrobacterium biovar 3) is the primary causative agent of crown gall disease of grapevine worldwide. We obtained and analyzed whole-genome sequences of diverse All. vitis strains to get insights into their diversification and taxonomy. RESULTS: Pairwise genome comparisons and phylogenomic analysis of various All. vitis strains clearly indicated that All. vitis is not a single species, but represents a species complex composed of several genomic species. Thus, we emended the description of All. vitis, which now refers to a restricted group of strains within the All. vitis species complex (i.e. All. vitis sensu stricto) and proposed a description of a novel species, All. ampelinum sp. nov. The type strain of All. vitis sensu stricto remains the current type strain of All. vitis, K309T. The type strain of All. ampelinum sp. nov. is S4T. We also identified sets of gene clusters specific to the All. vitis species complex, All. vitis sensu stricto and All. ampelinum, respectively, for which we predicted the biological function and infer the role in ecological diversification of these clades, including some we could experimentally validate. All. vitis species complex-specific genes confer tolerance to different stresses, including exposure to aromatic compounds. Similarly, All. vitis sensu stricto-specific genes confer the ability to degrade 4-hydroxyphenylacetate and a putative compound related to gentisic acid. All. ampelinum-specific genes have putative functions related to polyamine metabolism and nickel assimilation. Congruently with the genome-based classification, All. vitis sensu stricto and All. ampelinum were clearly delineated by MALDI-TOF MS analysis. Moreover, our genome-based analysis indicated that Allorhizobium is clearly separated from other genera of the family Rhizobiaceae. CONCLUSIONS: Comparative genomics and phylogenomic analysis provided novel insights into the diversification and taxonomy of Allorhizobium vitis species complex, supporting our redefinition of All. vitis sensu stricto and description of All. ampelinum. Our pan-genome analyses suggest that these species have differentiated ecologies, each relying on specialized nutrient consumption or toxic compound degradation to adapt to their respective niche.

Comparative Genomics of Novel Agrobacterium G3 Strains Isolated From the International Space Station and Description of Agrobacterium tomkonis sp. nov.

Strains of Agrobacterium genomospecies 3 (i.e., genomovar G3 of the Agrobacterium tumefaciens species complex) have been previously isolated from diverse environments, including in association with plant roots, with algae, as part of a lignocellulose degrading community, from a hospital environment, as a human opportunistic pathogen, or as reported in this study, from a surface within the International Space Station. Polyphasic taxonomic methods revealed the relationship of Agrobacterium G3 strains to other Agrobacterium spp., which supports the description of a novel species. The G3 strains tested (n = 9) were phenotypically distinguishable among the strains from other genomospecies of the genus Agrobacterium. Phylogenetic analyses of the 16S rRNA gene, gyrB gene, multi-locus sequence analysis, and 1,089-gene core-genome gene concatenate concur that tested G3 strains belong to the Agrobacterium genus and they form a clade distinct from other validly described Agrobacterium species. The distinctiveness of this clade was confirmed by average nucleotide identity (ANI) and in silico digital DNA-DNA hybridization (dDDH) comparisons between the G3 tested strains and all known Agrobacterium species type strains, since obtained values were considerably below the 95% (ANI) and 70% (dDDH) thresholds used for the species delineation. According to the core-genome phylogeny and ANI comparisons, the closest relatives of G3 strains were Agrobacterium sp. strains UGM030330-04 and K599, members of a novel genomospecies we propose to call genomovar G21. Using this polyphasic approach, we characterized the phenotypic and genotypic synapomorphies of Agrobacterium G3, showing it is a bona fide bacterial species, well separated from previously named Agrobacterium species or other recognized genomic species. We thus propose the name Agrobacterium tomkonis for this species previously referred to as Agrobacterium genomospecies 3. The type strain of A. tomkonis is IIF1SW-B1T (= LMG 32164 = NRRL B-65602). Comparative genomic analysis show A. tomkonis strains have species-specific genes associated with secretion of secondary metabolites, including an exopolysaccharide and putative adhesins and resistance to copper. A. tomkonis specific gene functions notably relate to surface adhesion and could be involved to colonize nutrient-poor and harsh habitats. The A. tomkonis strains from the ISS showed presence of a 40-kbp plasmid and several other potential mobile genetic elements detected that could also be part of conjugative elements or integrated prophages.

Different evolutionary trends form the twilight zone of the bacterial pan-genome.

The pan-genome is defined as the combined set of all genes in the gene pool of a species. Pan-genome analyses have been very useful in helping to understand different evolutionary dynamics of bacterial species: an open pan-genome often indicates a free-living lifestyle with metabolic versatility, while closed pan-genomes are linked to host-restricted, ecologically specialized bacteria. A detailed understanding of the species pan-genome has also been instrumental in tracking the phylodynamics of emerging drug resistance mechanisms and drug-resistant pathogens. However, current approaches to analyse a species' pan-genome do not take the species population structure into account, nor do they account for the uneven sampling of different lineages, as is commonplace due to over-sampling of clinically relevant representatives. Here we present the application of a population structure-aware approach for classifying genes in a pan-genome based on within-species distribution. We demonstrate our approach on a collection of 7500 Escherichia coli genomes, one of the most-studied bacterial species and used as a model for an open pan-genome. We reveal clearly distinct groups of genes, clustered by different underlying evolutionary dynamics, and provide a more biologically informed and accurate description of the species' pan-genome.

Phylogenomics reveals the basis of adaptation of Pseudorhizobium species to extreme environments and supports a taxonomic revision of the genus.

The family Rhizobiaceae includes many genera of soil bacteria, often isolated for their association with plants. Herein, we investigate the genomic diversity of a group of Rhizobium species and unclassified strains isolated from atypical environments, including seawater, rock matrix or polluted soil. Based on whole-genome similarity and core genome phylogeny, we show that this group corresponds to the genus Pseudorhizobium. We thus reclassify Rhizobium halotolerans, R. marinum, R. flavum and R. endolithicum as P. halotolerans sp. nov., P. marinum comb. nov., P. flavum comb. nov. and P. endolithicum comb. nov., respectively, and show that P. pelagicum is a synonym of P. marinum. We also delineate a new chemolithoautotroph species, P. banfieldiae sp. nov., whose type strain is NT-26T (=DSM 106348T=CFBP 8663T). This genome-based classification was supported by a chemotaxonomic comparison, with increasing taxonomic resolution provided by fatty acid, protein and metabolic profiles. In addition, we used a phylogenetic approach to infer scenarios of duplication, horizontal transfer and loss for all genes in the Pseudorhizobium pangenome. We thus identify the key functions associated with the diversification of each species and higher clades, shedding light on the mechanisms of adaptation to their respective ecological niches. Respiratory proteins acquired at the origin of Pseudorhizobium were combined with clade-specific genes to enable different strategies for detoxification and nutrition in harsh, nutrient-poor environments.

Phylogenetic distribution and evolutionary dynamics of nod and T3SS genes in the genus Bradyrhizobium.

Bradyrhizobium are abundant soil bacteria and the major symbiont of legumes. The recent availability of Bradyrhizobium genome sequences provides a large source of information for analysis of symbiotic traits. In this study, we investigated the evolutionary dynamics of the nodulation genes (nod) and their relationship with the genes encoding type III secretion systems (T3SS) and their effectors among bradyrhizobia. Based on the comparative analysis of 146 Bradyrhizobium genome sequences, we identified six different types of T3SS gene clusters. The two predominant cluster types are designated RhcIa and RhcIb and both belong to the RhcI-T3SS family previously described in other rhizobia. They are found in 92/146 strains, most of them also containing nod genes. RhcIa and RhcIb gene clusters differ in the genes they carry: while the translocon-encoding gene nopX is systematically found in strains containing RhcIb, the nopE and nopH genes are specifically conserved in strains containing RhcIa, suggesting that these last two genes might functionally substitute nopX and play a role related to effector translocation. Phylogenetic analysis suggests that bradyrhizobia simultaneously gained nod and RhcI-T3SS gene clusters via horizontal transfer or subsequent vertical inheritance of a symbiotic island containing both. Sequence similarity searches for known Nop effector proteins in bradyrhizobial proteomes revealed the absence of a so-called core effectome, i.e. that no effector is conserved among all Bradyrhizobium strains. However, NopM and SUMO proteases were found to be the main effector families, being represented in the majority of the genus. This study indicates that bradyrhizobial T3SSs might play a more significant symbiotic role than previously thought and provides new candidates among T3SS structural proteins and effectors for future functional investigations.

Diversity and evolution of surface polysaccharide synthesis loci in Enterobacteriales.

Bacterial capsules and lipopolysaccharides are diverse surface polysaccharides (SPs) that serve as the frontline for interactions with the outside world. While SPs can evolve rapidly, their diversity and evolutionary dynamics across different taxonomic scales has not been investigated in detail. Here, we focused on the bacterial order Enterobacteriales (including the medically relevant Enterobacteriaceae), to carry out comparative genomics of two SP locus synthesis regions, cps and kps, using 27,334 genomes from 45 genera. We identified high-quality cps loci in 22 genera and kps in 11 genera, around 4% of which were detected in multiple species. We found SP loci to be highly dynamic genetic entities: their evolution was driven by high rates of horizontal gene transfer (HGT), both of whole loci and component genes, and relaxed purifying selection, yielding large repertoires of SP diversity. In spite of that, we found the presence of (near-)identical locus structures in distant taxonomic backgrounds that could not be explained by recent exchange, pointing to long-term selective preservation of locus structures in some populations. Our results reveal differences in evolutionary dynamics driving SP diversity within different bacterial species, with lineages of Escherichia coli, Enterobacter hormaechei and Klebsiella aerogenes most likely to share SP loci via recent exchange; and lineages of Salmonella enterica, Citrobacter sakazakii and Serratia marcescens most likely to share SP loci via other mechanisms such as long-term preservation. Overall, the evolution of SP loci in Enterobacteriales is driven by a range of evolutionary forces and their dynamics and relative importance varies between different species.

Whole genome sequencing of Herpes Simplex Virus 1 directly from human cerebrospinal fluid reveals selective constraints in neurotropic viruses.

Herpes Simplex Virus type 1 (HSV-1) chronically infects over 70 per cent of the global population. Clinical manifestations are largely restricted to recurrent epidermal vesicles. However, HSV-1 also leads to encephalitis, the infection of the brain parenchyma, with high associated rates of mortality and morbidity. In this study, we performed target enrichment followed by direct sequencing of HSV-1 genomes, using target enrichment methods on the cerebrospinal fluid (CSF) of clinical encephalitis patients and from skin swabs of epidermal vesicles on non-encephalopathic patients. Phylogenetic analysis revealed high inter-host diversity and little population structure. In contrast, samples from different lesions in the same patient clustered with similar patterns of allelic variants. Comparison of consensus genome sequences shows HSV-1 has been freely recombining, except for distinct islands of linkage disequilibrium (LD). This suggests functional constraints prevent recombination between certain genes, notably those encoding pairs of interacting proteins. Distinct LD patterns characterised subsets of viruses recovered from CSF and skin lesions, which may reflect different evolutionary constraints in different body compartments. Functions of genes under differential constraint related to immunity or tropism and provide new hypotheses on tissue-specific mechanisms of viral infection and latency.

Two Novel Genomospecies in the Agrobacterium tumefaciens Species Complex Associated with Rose Crown Gall.

In this study, we explored the pathogenicity and phylogenetic position of Agrobacterium spp. strains isolated from crown gall tissues on annual, perennial, and ornamental plants in Iran. Of the 43 strains studied, 10 strains were identified as Allorhizobium vitis (formerly Agrobacterium vitis) using the species-specific primer pair PGF/PGR. Thirty-three remaining strains were studied using multilocus sequence analysis of four housekeeping genes (i.e., atpD, gyrB, recA, and rpoB), from which seven strains were identified as A. larrymoorei and one strain was identified as A. rubi (Rer); the remaining 25 strains were scattered within the A. tumefaciens species complex. Two strains were identified as genomospecies 1 (G1), seven strains were identified as A. radiobacter (G4), seven strains were identified as A. deltaense (G7), two strains were identified as A. nepotum (G14), and one strain was identified as "A. viscosum" (G15). The strains Rnr, Rnw, and Rew as well as the two strains OT33 and R13 all isolated from rose and the strain Ap1 isolated from apple were clustered in three atypical clades within the A. tumefaciens species complex. All but eight strains (i.e., Nec10, Ph38, Ph49, fic9, Fic72, R13, OT33, and Ap1) were pathogenic on tomato and sunflower seedlings in greenhouse conditions, whereas all but three strains (i.e., fic9, Fic72, and OT33) showed tumorigenicity on carrot root discs. The phylogenetic analysis and nucleotide diversity statistics suggested the existence of two novel genomospecies within the A. tumefaciens species complex, which we named "G19" and "G20." Hence, we propose the strains Rew, Rnw, and Rnr as the members of "G19" and the strains R13 and OT33 as the members of G20, whereas the phylogenetic status of the atypical strain Ap1 remains undetermined.

Coevolution of Sites under Immune Selection Shapes Epstein-Barr Virus Population Structure.

Epstein-Barr virus (EBV) is one of the most common viral infections in humans and persists within its host for life. EBV therefore represents an extremely successful virus that has evolved complex strategies to evade the host's innate and adaptive immune response during both initial and persistent stages of infection. Here, we conducted a comparative genomics analysis on 223 whole genome sequences of worldwide EBV strains. We recover extensive genome-wide linkage disequilibrium (LD) despite pervasive genetic recombination. This pattern is explained by the global EBV population being subdivided into three main subpopulations, one primarily found in East Asia, one in Southeast Asia and Oceania, and the third including most of the other globally distributed genomes we analyzed. Additionally, sites in LD were overrepresented in immunogenic genes. Taken together, our results suggest that host immune selection and local adaptation to different human host populations has shaped the genome-wide patterns of genetic diversity in EBV.

Minimal standards for the description of new genera and species of rhizobia and agrobacteria.

Herein the members of the Subcommittee on Taxonomy of Rhizobia and Agrobacteria of the International Committee on Systematics of Prokaryotes review recent developments in rhizobial and agrobacterial taxonomy and propose updated minimal standards for the description of new species (and genera) in these groups. The essential requirements (minimal standards) for description of a new species are (1) a genome sequence of at least the proposed type strain and (2) evidence for differentiation from other species based on genome sequence comparisons. It is also recommended that (3) genetic variation within the species is documented with sequence data from several clearly different strains and (4) phenotypic features are described, and their variation documented with data from a relevant set of representative strains. Furthermore, it is encouraged that information is provided on (5) nodulation or pathogenicity phenotypes, as appropriate, with relevant gene sequences. These guidelines supplement the current rules of general bacterial taxonomy, which require (6) a name that conforms to the International Code of Nomenclature of Prokaryotes, (7) validation of the name by publication either directly in the International Journal of Systematic and Evolutionary Microbiology or in a validation list when published elsewhere, and (8) deposition of the type strain in two international culture collections in separate countries.

Oral microbiomes from hunter-gatherers and traditional farmers reveal shifts in commensal balance and pathogen load linked to diet.

Maladaptation to modern diets has been implicated in several chronic disorders. Given the higher prevalence of disease such as dental caries and chronic gum diseases in industrialized societies, we sought to investigate the impact of different subsistence strategies on oral health and physiology, as documented by the oral microbiome. To control for confounding variables such as environment and host genetics, we sampled saliva from three pairs of populations of hunter-gatherers and traditional farmers living in close proximity in the Philippines. Deep shotgun sequencing of salivary DNA generated high-coverage microbiomes along with human genomes. Comparing these microbiomes with publicly available data from individuals living on a Western diet revealed that abundance ratios of core species were significantly correlated with subsistence strategy, with hunter-gatherers and Westerners occupying either end of a gradient of Neisseria against Haemophilus, and traditional farmers falling in between. Species found preferentially in hunter-gatherers included microbes often considered as oral pathogens, despite their hosts' apparent good oral health. Discriminant analysis of gene functions revealed vitamin B5 autotrophy and urease-mediated pH regulation as candidate adaptations of the microbiome to the hunter-gatherer and Western diets, respectively. These results suggest that major transitions in diet selected for different communities of commensals and likely played a role in the emergence of modern oral pathogens.

Ancestral Genome Estimation Reveals the History of Ecological Diversification in Agrobacterium.

Horizontal gene transfer (HGT) is considered as a major source of innovation in bacteria, and as such is expected to drive adaptation to new ecological niches. However, among the many genes acquired through HGT along the diversification history of genomes, only a fraction may have actively contributed to sustained ecological adaptation. We used a phylogenetic approach accounting for the transfer of genes (or groups of genes) to estimate the history of genomes in Agrobacterium biovar 1, a diverse group of soil and plant-dwelling bacterial species. We identified clade-specific blocks of cotransferred genes encoding coherent biochemical pathways that may have contributed to the evolutionary success of key Agrobacterium clades. This pattern of gene coevolution rejects a neutral model of transfer, in which neighboring genes would be transferred independently of their function and rather suggests purifying selection on collectively coded acquired pathways. The acquisition of these synapomorphic blocks of cofunctioning genes probably drove the ecological diversification of Agrobacterium and defined features of ancestral ecological niches, which consistently hint at a strong selective role of host plant rhizospheres.

Islands of linkage in an ocean of pervasive recombination reveals two-speed evolution of human cytomegalovirus genomes.

Human cytomegalovirus (HCMV) infects most of the population worldwide, persisting throughout the host's life in a latent state with periodic episodes of reactivation. While typically asymptomatic, HCMV can cause fatal disease among congenitally infected infants and immunocompromised patients. These clinical issues are compounded by the emergence of antiviral resistance and the absence of an effective vaccine, the development of which is likely complicated by the numerous immune evasins encoded by HCMV to counter the host's adaptive immune responses, a feature that facilitates frequent super-infections. Understanding the evolutionary dynamics of HCMV is essential for the development of effective new drugs and vaccines. By comparing viral genomes from uncultivated or low-passaged clinical samples of diverse origins, we observe evidence of frequent homologous recombination events, both recent and ancient, and no structure of HCMV genetic diversity at the whole-genome scale. Analysis of individual gene-scale loci reveals a striking dichotomy: while most of the genome is highly conserved, recombines essentially freely and has evolved under purifying selection, 21 genes display extreme diversity, structured into distinct genotypes that do not recombine with each other. Most of these hyper-variable genes encode glycoproteins involved in cell entry or escape of host immunity. Evidence that half of them have diverged through episodes of intense positive selection suggests that rapid evolution of hyper-variable loci is likely driven by interactions with host immunity. It appears that this process is enabled by recombination unlinking hyper-variable loci from strongly constrained neighboring sites. It is conceivable that viral mechanisms facilitating super-infection have evolved to promote recombination between diverged genotypes, allowing the virus to continuously diversify at key loci to escape immune detection, while maintaining a genome optimally adapted to its asymptomatic infectious lifecycle.

Ecological speciation in bacteria: reverse ecology approaches reveal the adaptive part of bacterial cladogenesis.

In this review, we synthesise current models and recent comparative genomic studies describing how bacterial species may emerge through adaptation to a new ecological niche and maintain themselves in the same niche over long time periods. We notably consider the impact of genetic exchange with phylogenetically close relatives living in sympatry and how this leads to the heterogeneous evolution of different genes within the bacterial genome. This heterogeneity provides landmarks to recognise genes that determine adaptation to the ecological niche, and we present reverse ecology strategies to unravel ecological properties of bacterial populations.

Inferring horizontal gene transfer.

Horizontal or Lateral Gene Transfer (HGT or LGT) is the transmission of portions of genomic DNA between organisms through a process decoupled from vertical inheritance. In the presence of HGT events, different fragments of the genome are the result of different evolutionary histories. This can therefore complicate the investigations of evolutionary relatedness of lineages and species. Also, as HGT can bring into genomes radically different genotypes from distant lineages, or even new genes bearing new functions, it is a major source of phenotypic innovation and a mechanism of niche adaptation. For example, of particular relevance to human health is the lateral transfer of antibiotic resistance and pathogenicity determinants, leading to the emergence of pathogenic lineages. Computational identification of HGT events relies upon the investigation of sequence composition or evolutionary history of genes. Sequence composition-based ("parametric") methods search for deviations from the genomic average, whereas evolutionary history-based ("phylogenetic") approaches identify genes whose evolutionary history significantly differs from that of the host species. The evaluation and benchmarking of HGT inference methods typically rely upon simulated genomes, for which the true history is known. On real data, different methods tend to infer different HGT events, and as a result it can be difficult to ascertain all but simple and clear-cut HGT events.